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ABSTRACT 





System health management based on condition monitoring is of great value and significance for improving 
system performance and reliability. As more data from condition monitoring is obtained, the efficiency of system 
health monitoring improves. However, increasing the amount of data by installing new sensors on system 
components is restricted by a variety of factors such as budget, weight and space allowance. Therefore, sensor 
selection in the most optimal manner in order to achieve maximum valuable information from system is an 
important issue. In this study, a novel concept is introduced that can be used in sensor combination optimization 
problems. This concept employs failure modes and effect analysis method in signal information value quantifi- 
cation procedure. 

To verify the effectiveness of the proposed methodology, an optimization tool is developed based on the 
proposed methodology for an air handling unit as a case study. This tool offers sensor combination optimization 
and visualization of key failure mode parameters. Results show that the optimal sensor combination lead to a 
maximum information on system health considering budget limits. In addition, the sensitivity analysis on 
severity factor of component failures is performed. Results indicate that, by increasing the severity factor of a 
component failure, the optimization model select sensors that monitor more information on that failure mode. 
According to the results obtained, the proposed methodology can be used in maximizing the performance of 
health monitoring and management of complex systems with numerous applications including aerospace, oil and 


gas industry, and medical technology. 





1. Introduction 


Currently, complex energy systems may use a variety of primary 
energy sources and different equipment to supply building energy de- 
mand. This goal is achieved by various operating modes of equipment. It 
is crucial to determine the best operating mode of an equipment 
involved in such a way that the energy consumption of the system can be 
minimized [1]. Many studies have focused on the relation between 
system energy consumption and operating conditions [2,3]. The energy 
consumption as a function of operating conditions can be derived using 
principle-based and/or data-driven methods [4]. 

In addition to the operating conditions, malfunction of equipment 
can cause abnormal energy consumption. Root cause analysis of 
equipment malfunctions in order to reduce energy consumption is a 
critical issue [5]. Fault detection and isolation in complex systems has 
been accompanied by development of advanced algorithms, applica- 
tions and services. A review on computational intelligence in fault 
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detection and energy management can be found in literature [6]. 

In order to diagnose and isolate a fault in a system, a sensor network 
(SN) should be implemented on the system. SNs can detect useful in- 
formation about system operating status [7]. In recent years, there has 
been an increasing interest in the use of SNs for system health moni- 
toring and management [8]. SNs consist of a network of sensor devices, 
where each device can sense the target operational and environmental 
conditions to achieve the goal of delivering valuable information to the 
health monitoring system. Various factors should be considered when 
deploying a SN, including placement and combination selection [9]. 

For the placement, the location of a specific sensor type is optimized 
in order to meet a certain criterion [10]. Sensor placement optimization 
is widely applied in water distribution networks. The optimization 
models determine the minimum number of sensors and their optimal 
placement in order to satisfy a pre-specified criterion for the water 
network such as reliability of contaminant detection [11]. Recently, 
there has been growing interest in sensor placement optimization in 
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energy field as well. In Yoganathan et al. [12] study, the location and 
number of three sensor types were optimized separately. A data-driven 
approach was used for the optimization procedure, and the optimal 
sensor network provided a holistic picture of the indoor environment. 

There are several studies on sensor placement optimization [13-15]. 
However, research on sensor combination optimization is limited. In 
sensor combination optimization, the goal is to find the optimal types of 
sensors (e.g. thermometer, flowmeter, pressure meter) in a complex 
system with multiple target operational and environmental parameters. 

As is clear, it is not possible to have sensors on all target operational 
and environmental parameters due to the lack of budget or space/ 
weight limits - i.e., some engineering systems have weight and space 
restrictions in the design phase. For instance, weight and size limitations 
on the spacecrafts are severe. As a result, the number of sensors that 
could be installed are limited. Another example is heating, ventilation, 
and air-conditioning (HVAC) systems with limited available space. Most 
of the residential HVAC packaged units have a limited available space to 
install additional sensors. Therefore, its critical to use this space effec- 
tively, and install a sensor network that provide maximum information 
about system health and performance. 

As mentioned, optimal sensor combination is a critical design issue, 
as it directly affects the performance of the system health monitoring 
system. However, few studies have focused on this topic. In Parhizkar 
et al. [9], the comparison between leakage and corrosion sensors in an 
oil pipeline was performed. The optimal Pareto front of optimal sensor 
combinations as a function of budget limit was presented and discussed 
in this study. 

In optimization of sensor combination problems, the value of signals 
from different sensor types should be quantified in a manner that results 
be comparable. The value of a signal depends on its ability to detect an 
event or fault in a system. The information value concept is introduced 
in this study in order to quantify the signal values of a SN. 

Information value is defined as the value of information that is 
conveyed by a sensor network to detect system failure modes and health 
status. This value shows how efficient a SN can detect the status of all 
events that are probable to happen in a system. In order to compare 
different sensor combinations, the information value should be quanti- 
fied. There can be several methods for quantifying information value. 
One of the most efficient methods of qualification can be failure modes 
and effect analysis (FMEA) that is investigated in this study. 

The sensor combination optimization (SCO) based on FMEA method 
has the potential for being used to optimize sensor combination of any 
complex energy system. 

Across the United States, the average annual energy intensity for an 
office building is 906 kW per square meter and the average cost is $189 
per square meter [16]. The energy consumption on space HVAC repre- 
sents about 55% of the total energy use in a typical office building [17]. 
The air handling unit (AHU) of an HVAC system impacts the overall 
performance of the system. The goal of this study is to minimize the total 
energy consumption of an AHU system while detecting system faults 
with highest accuracy at different budget limits for monitoring system. 
The proposed method is applied on an AHU system and the performance 
of the method is discussed. The main contributions of this research can 
be summarized as follows: 


1. This study illustrates the application of fault detection for energy 
consumption reduction by detecting upcoming faults automatically 
and eliminating unexpected failure modes. 

2. As mentioned, research on sensor combination optimization is not 
well developed. In this study a new method for sensor combination 
optimization is proposed. In sensor combination optimization, type 
of sensors is optimized in addition to sensors location. The proposed 
method is based on FMEA, and it is applied to a case study. 

3. A graphical sensor optimization and fault detection software is 
developed for design and operation phases of AHU systems. The 
software is able to optimize sensor combination at the design phase. 
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In addition, system failure modes and abnormal energy consumption 
could be detected at the operating phase of the AHU system. In the 
following sections, the main concept of the proposed framework and 
the methodology of the optimization are described. Also, the appli- 
cation of the proposed method on an AHU system and the results are 
discussed and the optimal sensor combination and its fault detection 
performance are presented. 


2. Main concept of proposed framework 


Health monitoring of complex energy systems require multiple sen- 
sors to extract required information from the sensed internal conditions 
and environment of the system and its components [18]. A critical de- 
cision, particularly in the context of complex energy systems, is the type 
of the sensors given a set of technical and non-technical constraints. The 
proposed methodology is designed to answer important questions such 
as how to infer system health status based on limited number of moni- 
toring parameters at certain sub systems. The concept of the proposed 
methodology is presented in Fig. 1. 

As can be seen in Fig. 1, ina limited budget and/or weight allowance, 
different sensor combinations can be selected. However, each sensor 
combination detects different parameters and, as a result, the informa- 
tion about system health status is variable. In this paper, a method is 
proposed that can find the optimal sensor combination which leads to a 
maximum system health information in a restricted condition (budget 
and weight limits). 


3. Sensor combination optimization methodology 


This paper provides a FMEA-based SCO methodology for health 
monitoring of energy systems. The approach uses the physical models of 
the system, and FMEA techniques along with the constraints. FMEA is a 
systematic methodology to quantify the risk and failure of engineered 
systems. A FMEA combines the failure modes, the causes of failures, the 
effects of failures, and corrective actions [19]. 

The flow diagram of the proposed methodology is presented in Fig. 2. 
By casting the FMEA in a probabilistic framework through a measure of 
the likelihood of occurrence of each failure mode and a measure of the 
consequence of each failure mode severity of the effect, the risk asso- 
ciated with each failure mode can be evaluated. 

In addition, each sensor combination has a failure detection perfor- 
mance that can be calculated through failure detection module, and the 
sensor combinations are ranked. The ranking is based on the Modified 
Risk Priority Number (MRPN), which is the product of the measure of 
likelihood, severity, and detection performance. A MRPN is calculated 
for each sensor combination, and subsequently by considering available 
constraints, the optimal configuration is determined through optimiza- 
tion module. 


3.1. Expert rule/failure tree generation 


In this step, the possible failure modes are identified. The list should 
be exhaustive and include all possible failure modes. Those perceived as 
not likely to occur or with negligible consequence can be eliminated 
from further consideration at this stage [20]. For those remaining failure 
modes, the FMEA process will quantify the risk associated with each 
failure mode using the MRPN. This process will allow for elimination of 
those failure modes that are considered insignificant. 


3.2. Failure mode probability 


In this step, the probability of occurrence of each failure mode and 
the associated likelihood measure scale is developed. The scale of like- 
lihood can be developed based on values or quantitative measures such 
as statistics, as is often the case for qualitative assessment [21]. 
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Fig. 1. Sensor combination optimization concept. 
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Fig. 2. Flow diagram of the FMEA-based sensor combination optimization. 


3.3. Failure mode severity 


In this step, a scale for severity of the consequences associated with 
each failure mode is developed. The consequences are based on quali- 
tative (e.g. level of erosion) or quantitative (e.g. cost to repair) metrics 
[21]. 


3.4. Failure mode detection 


The detection performance of a sensor combination is calculated in 
this step. In this study, the performance of a SN in detecting events status 


is calculated using a detection matrix. The detection matrix presents all 
events or failure modes in rows and sensor types in column. For instance, 
for an experts rule tree illustrated in Fig. 3, the detection matrix has 4 
rows and 3 columns. 

The row elements of the detection matrix can be equal to "a" or 1. The 
required sensors to detect an event of a row is presented as"a" and the 
rest is equal tol. The detection matrix of the event of Fig. 3 is equal to: 
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Fig. 3. Experts rule tree sample. 


T m P 
ar am 1 
ar dm 1 (1) 


IM = 


The presented matrix is a general form of an information matrix. The 
information matrix of a system is equal to the number of sensor com- 
binations that can be calculated as Eq. (2). 

C(n,k 4 (2) 
n <= —ÀÁáÁI 
(n, ) k\(n — k)! 

For instance, the presented example in Fig. 3, has 7 information 

matrices. 


C(3,1)--C,2)--C(3,3) = 34-3341 = 7 (3) 


Each sensor combination has a unique information matrix that the 
“a” values are equal tol or 0.5, depending on the presence of a sensor in 
the combination. If sensor type j is included in the combination, then all 
the jth column with “a” value are equal to 1, and if the sensor type j is not 
in the sensor combination, all “a” values in the jth column are equal to 
0.5. Therest of the matrix elements with value 1, remain constant. In the 
presented example with thermometer and mass flowmeter, the infor- 
mation matrix is: 


T m P 

1 

1 (4) 
0.5 


1 
1 
IM = 1 
1 0.5 


1 
1 
1 
1 

The product of elements listed in the ith row, represents the infor- 
mation value of that sensor combination in detecting ith event/failure 
mode. The information values of the presented example with a ther- 
mometer and a mass flowmeter is as Eq. (5). As can be seen, each row 


has one element that is the product of all elements in the same row of IM 
matrix. 


1 
IV = 0.5 (5) 


0.5 


As the information value is the product of row elements of infor- 
mation matrix, it is unitless and the lower and upper bands are equal to 
0 and 1, respectively. 


3.5. Optimization model 


The objective of the optimization model is to maximize the overall 
information value of a sensor network. In this study, it is assumed that 
the overall information value is equal to the modified risk priority 
number (MRPN). 


Max MRPN (6) 


Where, MRPN is equal to the sum of the overall information value (OIV) 
arrays and is calculated as Eq. (7). 
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MRPN = X OIV, (7) 
i=1 
OIV is the overall information value, which is equal to the product of 
the information value matrix, probability, and severity arrays, and is 
calculated as Eq. (8). 


OIV=PxSxIV (8) 


P is the failure probability, S is the failure severity, and IV is the 
information value. Therefore, the objective function could be presented 
as Eq. (9). 


Max MRPN = >. OIV; = NOP, x S; x IV; (9) 
— > 


i=l 


l= 


The constraints of the optimization model could be weight, space, or 
budget limits of the SN, according to its specific application. For 
instance, in aerospace SNs have weight limit restrictions [22], or in oil 
industry budget limits are prevalent [23]. For n sensors, the weight and 
budget limits are as Eqs. 10 and 11, respectively. 


Som < W' (10) 
i=1 
Ys < C (11) 


i=l 
4. Application of the proposed method on an ahu system 
4.1. System description 


AHUs are used to modify and circulate air to satisfy comfort condi- 
tions in buildings. AHUs are mainly composed of two fan coils and a 
damper connected through a large case. AHU usually connects to a 
ductwork ventilation system that distributes the conditioned air through 
the building and returns it to the AHU [24]. Fig. 4 illustrates the sche- 
matic of an AHU and the parameters that can be detected by sensor 
networks. 

In addition to the presented parameters in Fig. 4, room and ambient 
air physical properties are required for system fault detection and energy 
consumption analysis. All these parameters can be detected by three 
sensor types including temperature sensor, relative humidity sensor, and 
flow meter sensor. All the required sensors for AHUs fault detection and 
energy consumption analysis are identified [25], and presented in 
Table 1. 

The case study selected for the analysis, in this study, is the AHU of 
an office building located in Tehran, Iran. The main characteristics of the 
studied AHU is presented in Table 2. More details about the system is 
provided in Ref. [26]. 


4.2. System model 


Sensor network can be any combinations of presented sensors in 
Table 1. As mentioned in the previous section, the first step is fault/ 
experts rule tree generation. 

For deriving the experts rule tree, the operation of the AHU should be 
identified, and, thermodynamic rules such energy and mass balance for 
all parts of AHU should be derived. The AHU has four operational modes 
including cooling mode with maximum inlet air flow from outside, 
cooling mode with minimum inlet air flow from outside, no heating and 
cooling mode, and heating mode. In this research, the first two modes 
are merged [27], and three main operating modes of cooling (mode 1), 
no cooling/heating (mode 2), and heating (mode 3) are considered in 
the analysis. 

In mode 1, cooling coil is open and heating coil is closed. The amount 
of inlet air flow from outside is variable to satisfy the comfort condition 
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Fig. 4. Detectable parameters by sensor network in an AHU system. 


Table 1 
Required sensors in an AHU system for fault detection and energy consumption 
analysis. 


Sensor number Symbol Sensor description 

1 To Outdoor thermometer 

2 T; Return air thermometer 

3 Tsup Supply air thermometer 

4 T Coils interconnection thermometer 
5 Taig Mixed air thermometer 

6 Tewa Inlet cold water thermometer 
7 Tewo Outlet cold water thermometer 
8 Thw,i Inlet hot water thermometer 

9 Thw,o Outlet hot water thermometer 
10 Vo Outlet air flowmeter 

11 Vr Return air flowmeter 

12 Vs Supply air flowmeter 


Hot water flowmeter 
Cold water flowmeter 


= me 
AC 
HE 
= š 


15 RH, Outdoor air relative humidity sensor 

16 RHgix Mixed air relative humidity sensor 

17 RH, Return air relative humidity sensor 
Table 2 


Characteristics of the studied AHU system. 


Characteristic Description 


Air flow range (m?/h) 1000-100000 


Basic options e G4 class filters 
e F5,F7,F9 class filters 
e Humidification, fan and attenuation section 


Insulation type Indoors (YMBS)/outdoors (YMBD) 


and decrease energy loss. In mode 2, heating and cooling coils are closed 
and the comfort condition are satisfied by the outdoor air. In mode 3, 
heating coil is open and cooling coil is closed. The inlet air flow from 
outside is minimum to decrease waste of energy. 

After recognizing the operational mode of the AHU, thermodynamics 
rules should be studied to infer failure modes in the system. The laws of 
thermodynamics for each operating mode are presented in Table 3, 
[28-30]. 

The experts rule tree of the AHU is derived using presented rules in 
Table 3, that is illustrated in Fig. 5. Tree based learning algorithms are 
considered to be one of the best and mostly used supervised learning 


Table 3 
Thermodynamics rules for failure mode detection of the AHU [28-30]. 
Operational mods Rules 
Cooling mode (mode 1) ITs« — Ts| < € 
To > Tg IT — Ts| < € 


|T rie m TMA idea | <E 
Vg ca — Vr] < € 





Im, = Myc_cal| < E 
No heating and cooling mode (mode 2) [Tse — Ts| « € 
Ts « To « Tg IT. — Ts| < € 


(Trie — Tua ideat| < € 





IVa ca m Val <E 
m, —hn,. cal «e 
Heating mode (mode 3) ITs« — Ts| < € 
To < Tet Ty = Ts| <E 


Tmix — TMA ideal | <E 





Myn — Myh_cal| < E 


Vp Vile 


-cal 





Vo. cat n Vol <E 


methods. Tree based methods empower predictive models with high 
accuracy, stability and ease of interpretation. Unlike linear models, they 
map non-linear relationships quite well. They are adaptable at solving 
any kind of problem at hand (classification or regression). An expert rule 
tree is a set of technical rules that uses a tree-like graph or model of rules 
and their possible consequences. It is one way to display an algorithm 
that only contains conditional control statements. 

An expert rule tree is a flowchart-like structure in which each in- 
ternal node represents a “test” on a variable (e.g. whether temperature is 
in a specific range or not), each branch represents the outcome of the 
test, and each leaf node represents a class label (system status after 
computing all variables). The paths from root to leaf in Figs. 5-7 
represent classification technical rules for modes 1, 2 and 3, respec- 
tively. As can be seen, expert rule trees cover critical control policies as 
well. 

According to Figs. 5-7, failure modes of the system is presented in 
Table 4. 

Failure mode probabilities are given in Ref. [31], and presented in 
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Fig. 5. The experts rule tree of the studied AHU for operating mode 1. 
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Fig. 6. The experts rule tree of the studied AHU for operating mode 2. 


Table 5. In addition, failure mode severity is assumed to be proportional 
to the average repair time of the related component, given in Ref. [31]. 
The severity value of the failure modes are presented in Table 5. 

In order to calculate overall information value, the product of 
probability and severity for all failure modes are calculated. The derived 
values are normalized to have values between 0 and 1. The normaliza- 
tion is performed using Eq. (12), and the normalized values are pre- 
sented in Table 6. 

x; — min(x) 


a max(x) — min(x) ied 


Next step is the failure mode detection, where information matrix is 
developed. 

As discussed, each row represents an event and the columns show the 
sensor types. The AHU information matrix has 27 rows, which presents 
27 events, and it has 9 columns for 9 sensor types. The required sensors 
for each event is presented as "a", and if the sensor is not required to 
detect the event, the value is equal to 1. 
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Fig. 7. The experts rule tree of the studied AHU for operating mode 3. 


Failure modes and sensor status of the AHU. 


Number 


1 
2 


NO OF Bp W 


10 
11 


12 


13 
14 


Table 5 


Failure mode 


Mode 1- Opening valve fault 
Mode 1- Leakage a 


Mode 1- Chiller fault 

Mode 1- Sensor error a 
Mode 1- Leakage b 

Mode 1- No fault 

Mode 1- Outdoor air damper 
fault 

Mode 1- Return air damper 
fault 

Mode 1- Sensor error b 
Mode 2- Opening valve fault 
Mode 2- Leakage a 


Mode 2- Chiller fault 


Mode 2- Leakage b 
Mode 2- Sensor error a 


Number 


15 
16 


17 
18 
19 
20 
21 


22 
23 
24 
25 


26 


27 


Failure rate and severity factor of failure modes. 


Failure mod 


Opening val 
Leakage 
Chiller fault 


Sensor error 
Damper error 


e Failure rate [31] 


ve fault 5E-06/h 
3.5E-06/h 
3E-06/h 


8E-06/h 
3E-06/h 


T -T 


x mix 








<E 


Yes 





Failure mode 


Mode 2- No fault 

Mode 2- Return air damper 
fault 

Mode 2- Sensor error b 
Mode 3- Sensor error a 
Mode 3- Opening valve fault 
Mode 3- Leakage a 

Mode 3- Boiler fault 


Mode 3- Leakage b 


Mode 3- Sensor error b 
Mode 3- No fault 

Mode 3- Outdoorair damper 
fault 

Mode 3- Return air damper 
fault 

Mode 3- Sensor error d 


Severity (h) [31] 


19 
12 
40 
10 
24 
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pipes 
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Normalized factors for considered failure modes. 


Failure mode number 


Opening valve fault 
Leakage 

Chiller fault 

Sensor error 
Damper error 


Normalized factor 


0.8 
0.35 
1 
0.69 
0.6 
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a l1 1 1 1 1 1] I a 
a 1 1 1 1 1 1 1] a 
a 1 1 a 1 1 1 1 1] 
a4 1 1 d4 1 1 1 dg 1 
a4 1 1 d4 1 1 1 dg 1 
1 1 1 1 1 1 a 1 1 
1 1 1 1. as 1 a 1 1 
1 1 1 1 ü5 ag Q3 1 1 
1 1 1 1 ds ag Q3 1 1 
a 1 1 1 1 1 1] 1] a 
a 1 1 1 1 1 1 1] a 
a 1 1 a 1 1 1 1 1] 
dı 1 1 d4 1 1 1 dg 1 
A= dı 1 1 d4 1 1 1 dg 1 (13) 
1 1 1 1 1 1I a 1 1 
1 1 1 1 1 a% a 1 1! 
1 1 1 1 1 a% a 1 1! 
1 L 21 1 Ltt 1 1 1 
aa 1 1 1 1 1 a 1 
aa 1 1 1 1 1 a 1 
d a a 1 1 1 1 1 1 
Qa, a Q3 1 1 1 1 1 ao 
Qa, a Q3 1 1 1 1 1 ao 
1 1 1 1 1 1] a 1 1! 
1 1 1 1] a5 1 a 1 1! 
1 1 1 1 ü5 agp Q3 1 1 
1 1 1 1 ü5 ag Q3 1 1 


The presented information matrix shows the general form of one 
sensor combination. It should be noted that the combinations of the 9 
sensor types is equal to 511 (Eq. (2)). As a result, there are 511 infor- 
mation matrices for 511 sensor combinations. In each sensor combina- 
tion, If sensor type j is included in the combination, then all the j th 
column with “a” value are equal to 1, and if sensor type j is not included 
in the sensor combination, then all “a” values in the j th column are 
equal to 0.5. As mentioned before, the product of the elements in i row, 
represents the information value of the sensor combination in detecting 
ith event. 

The final step is the optimization model. The objective function of 
the model is the maximization of the MRPN that is presented in Eq. (6). 
In addition, the model is subject to a budget constraint as Eq. (14). This 
constraint represents that the SN costs must be less than the available 
budget for investment. 


5% < C (14) 


i=l 


4.3. Solution algorithm 


The problem under study is a mixed-integer nonlinear programming 
(MINLP) problem. MINLP addresses a class of optimization problems 
with nonlinearities in the constraints and/or objective as well as integer 
and/or continuous variables. Multiple solution algorithms are proposed 
for MINLP problems. According to the performance of the solution al- 
gorithm, particle swarm optimization (PSO), which is the best suitable 
algorithm, is selected for the case under study. 

The PSO method is a population-based optimization algorithm 
inspired by social behavior of fish schooling and bird flocking. PSO al- 
gorithm shares many similarities with other evolutionary optimization 
algorithms such as Genetic Algorithms (GA). One of the main advantage 
of PSO algorithm over other evolutionary methods is that there are few 
parameters to be adjusted in PSO algorithm. The algorithm starts with a 
population of random solutions and searches for optimal point by 
updating generations. In other words, the PSO computational method 
optimizes a problem by iteratively improve a candidate solution with 
considering a given measure of quality. The candidate solution is 
influenced by the local best-known point and it is guided toward the 
best-known positions in the search-space, which are updated as better 
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positions are found by other particles. This process is expected to move 
the candidate solution toward the optimal solution [32]. 

The PSO algorithm is suitable for nonlinear programming with 
continuous decision variables. However, different algorithms are pro- 
posed to modify PSO for handling mixed integer nonlinear program- 
ming. One of the proposed modifications is to adapt random search 
methods to the integer case simply by adding a round off instruction to 
choose the nearest integer of a given real value, which is utilized in this 
study. 


4.4. Graphical user interface development 


The SCO model for AHU system is visualized in Matlab software. To 
ensure that the tool would be of maximum utility to operators, design 
feedback was obtained iteratively, throughout development. Fig. 8 
shows the landing page of the tool. In the plot, the total available budget 
should be given as an input. Next, by selecting the optimum combina- 
tion, the optimal number of each sensor type would be determined. For 
instance, by selecting $70 budget limit, the optimal sensor combination 
would be thermometers on outdoor, return air, supply air, outlet hot 
water, mixed air, inlet hot water, outlet cold water, inlet cold water, 
coils interconnection, and flowmeters on supply air, hot water, and cold- 
water streams. 

In the presented tool, the severity of each event and failure mode is 
pre-defined based on the historical data of AHUs in office buildings. 
However, this value can be dependent on the application of the AHU as 
well. For instance, AHU leakage in a manufacturing process has a 
different severity than it has in a residential building. 

In order to consider different severity factors, the tool is modified. 
The modified tool, in addition to the budget limit, receives the scaled 
severity value of failure modes as an input. The modified tool interface is 
presented in Fig. 9. The plot in the loading page shows the scaled in- 
formation value of failure modes for the optimal SN. Number ten in- 
dicates that the event/failure mode can be fully detected by the optimal 
sensor combination and zero value represents no information is obtained 
about the related event/failure mode. 


5. Results 
5.1. Model performance evaluation 


In this section, a sensitivity analysis to investigate the performance of 
the model is performed. As discussed, failure modes have different 
severity factors depending on the system application. This fact has an 
effect on the optimal sensor combination in system health monitoring. It 
is desirable to detect failure modes with higher severity factor. There- 
fore, as the severity factor of a failure mode increases, the optimization 
model should select sensors for detection of that failure mode. Table 7 
presents the optimal sensor combinations for four different scenarios. 
Scenarios 1 consider chiller fault severity factor equal to one, which 
means the most important component that its health status should be 
monitored in the system is chiller. In scenario 2, the severity factor of 
chiller fault is zero. In this scenario chiller function is not important in 
the system, that can be a reasonable scenario in cold seasons. Finally, 
scenario 3 and 4, consider damper fault severity equal to O and 1, 
respectively. 

Fig. 10 presents the information value of failure modes in scenario 1 
and 2. In scenario 1, optimal sensor combination monitor failure modes 
of 3, 12, 18, and 21, perfectly. According to Table 4, these failure modes 
are chiller fault, chiller fault, sensor error and boiler fault, respectively. 
It means that the maximum information on chiller function is monitored 
with considering the chiller fault severity factor equal to one. As severity 
factor of chiller fault decreases, optimal sensor combination is selected 
based on the monitoring more information about the rest of the system. 
As a result, the information value of multiple failure modes except 
chiller faults are high in this case, that is shown in Fig. 10 (b). 
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Fig. 8. Screenshot of the landing page of the SCOtool. 
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Fig. 9. Screenshot of the landing page of the modified SCO tool. 


The information value of failure modes for scenario 3 and 4 is pre- 
sented in Fig. 11. As can be seen, in damper fault severity factor equal to 
zero, system monitor failure modes of chiller fault completely. That is 


because of the severity factor of chiller fault, which is equal to one, 
according to Table 6. By increasing damper fault severity factor to one, 
system monitor more information about dampers health status. Fig. 11 
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11, higher severity factor of a component results in an optimal SN that 
monitor more information about that component. 


5.2. Information value of sensor combinations 


As mentioned before, the value of information of different sensor 
combinations, depends on the failure modes probabilities, severities, SN 
detection performance. Fig. 12 presents the information value of all 
possible combinations of sensors in the studied AHU versus SN cost. 

As can be seen in Fig. 12, there is an upward trend in information 
value. As the available budget increases, more sensors can be purchased, 
and as a result, more valuable information can be gathered. In addition, 
Fig. 12 shows that all sensor combinations are districted in six separated 
areas. That is because of the price ratio of sensors. As presented in 
Table 1, there are three type of sensor in the system including ther- 
mometer, flowmeter, and humidity sensor. The price of flowmeter is 
much higher that the two other sensors. As a flowmeter is added to the 
sensor combination, a spike in the SN cost will occur that is shown in the 
figure. 


Table 7 
Optimal sensor combination for different severity factors of chiller and damper fault. 

No 4 5 6 
1 Chiller fault severity 1 1 1 1 0 0 1 
2 Chiller fault severity O 1 1 1 1 1 0 
3 Damper fault severity 0 1 1 1 0 0 1 
4 1 1 1 0 0 1 
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(b) Scenario 2: Chiller fault severity=0 


Fig. 10. Sensitivity analysis of severity factor of chiller fault on optimal sensor 
combination. 
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(a) Scenario 3: Damper fault severity= 0 
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(b) Scenario 4: Damper fault severity- 1 


Fig. 11. Sensitivity analysis of severity factor of damper fault on optimal sensor 
combination. 


(b) shows that failure modes of 7, 8, 16, 18, and 25 are monitored 
completely. These failure modes are outdoor air damper fault, return air 
damper fault, return air damper fault, sensor error and outdoor air 
damper fault, respectively. In other word, as presented in Figs. 10 and 
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5.3. Optimal sensor combinations 


To obtain the optimal sensor combination, budget constraint is 
applied to the model. Fig. 13 depicts a budget constraint of $140. In this 
case, a sensor combination with maximum information value in the 
feasible area (shaded blue) is selected as the optimal solution. 

In other energy systems, that sensor costs are almost the same, the 
gap between sensor combination is much narrow, and as a results, the 
combination optimization is more complex. For instance, Fig. 14 shows 
the sensor combinations as a function of SN cost for a chiller system. 
Sensor types in a chiller system include flowmeter, thermometer, and 
pressure sensors, where the costs of thermometer and pressure sensor 
are almost in the same range. Therefore, sensor combinations are spread 
widely, and as a result, the optimal combination trend are not always 
ascending as a function of cost. As illustrated in Fig. 14, the optimal 
combination, with $140 budget limit, costs around $120, and more 
expensive combinations ($120 to $140) do not results in more valuable 
information. This fact shows the importance of optimal budget alloca- 
tion in sensor selection in order to get most efficient information about 
system health status. 

Table 8 presents the optimal combination of sensors at different 
budget limits. As is shown, at higher budget limit, number of sensors are 
increased. 





Overall information value (OIV) 
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Fig. 12. Information value as a function of cost for all sensor combinations in 
the studied AHU. 
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Fig. 13. Information value versus cost with considering a budget constraint for 
AHU system. 
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Fig. 14. Information value versus cost with considering a budget constraint for 
chiller system. 


5.4. Fault detection performance 


In this section, fault detection performance of the optimal sensor 
combinations is presented. Fig. 15 illustrates the gained information 
from SN about all system failure modes. The gained information is 
presented based on a scaled information value index. Higher index 
means more information is available about the system. It can be inferred 
from Fig. 15, by a minimum available budget, failure modes of 3, 12, 18, 
and 21 can be fully detected. According to Tables 3 and 4, these failure 
modes are all of the failure modes that require a flowmeter and multiple 
thermometer for detection. As the budget increases, the number of flow 
meters can be increased, and as a result more failure modes can be fully 
detected. Fig. 15(f) presents a case with a budget limit of $300. As the 
combination of all the sensors in this system is cheaper than $300, this 
constraint doesn't change the feasible region of optimization. As a 
result, the optimization model select all the available sensors as the 
optimal sensor combination, and this combination fully detects all fail- 
ure modes as presented in Fig. 15(f). 


5.5. Abnormal energy consumption detection 


This study described AHU failure modes for detecting abnormal 
energy consumption in AHU system. In fact, the standard required 


Table 8 
Optimal sensor combinations with considering different budget limits. 


Budget limit Optimal combination 


1 2 3 4 5 6 7 8 
< $50 1 1 1 1 1 1 1 1 
< $100 1 1 1 1 0 1 1 0 
« $150 1 1 1 1 1 1 1 0 
« $200 1 1 1 1 1 1 1 1 
« $250 1 1 1 1 0 1 1 1 
« $300 1 1 1 1 1 1 1 1 
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energy for system operation at the sensed operating and environmental 
conditions is compared with the actual energy consumption of the sys- 
tem. A acceptable deviation range should be defined. As signals show, 
actual energy consumption exceeds the acceptable deviation range from 
the standard value, the fault detection and isolation module would be 
activated to search for the main cause of such abnormalities, detect the 
event/fault, and isolate it as quickly as possible in order to reduce sys- 
tem abnormal energy consumption. Fig. 16 shows the deviation of actual 
energy consumption from calculated standard value for 100 random 
signals from sensors. The fault area shows the states that have an 
abnormal energy consumption. According to the developed experts rule 
tree and fault detection module, the possible failure modes that have 
caused abnormal energy consumption can be detected. As illustrated in 
Fig. 16, by decreasing the deviation range, more events are reported as 
possible predicted failure modes. In order to the system consume energy 
in a range closer to its standard value, more events are detected as the 
potential failure modes. In other words, energy consumption of the 
system can be reduced more efficiently by detecting and isolating more 
potential failure modes in the system. However, as the deviation range 
increases to 2596, system allowable energy consumption range increase, 
and most of the states operate in an acceptable energy consumption 
range. 

The detected faulty state is analyzed in the developed experts rule 
tree, and root causes of the possible failure modes that lead to energy 
loss is detected. For instance, the blue dot in Fig. 16 is a signal of [To, 
RHO, Tr, Tsup, Tset, Tx, Tmix, Twic, Twoc, Mwc, Twih, Twoh, Mwh, Vs, 
Vo, Vomin, Vr, RHmix, RHR]= [34, 16, 29, 35.5, 18, 21, 32.5, 31.3, 
34.6, 35.9, 0.1, 25.8, 25.5, 0.1, 0.7, 0.8, 0.4, 0.5, 0.6, 0.4] that is 
detected for the cause of abnormal energy consumption. By analyzing 
the signal through the proposed experts rule tree, as presented in Fig. 17, 
chiller fault is detected. 


6. Conclusion and recommendation 


Sensor network in health monitoring system should be optimized due 
to budget, space, and weight limits in complex systems. In this study, 
sensor combination optimization concept is introduced. This concept 
has the advantage over previous approaches in that it optimizes sensor 
types, in addition to sensor locations. This concept optimizes sensor type 
and location according to the gained information about system failure 
modes. The gained information value is calculated based on the failure 
modes and effect analysis, and optimal sensor combination maximizes 
the performance of the health monitoring system. The failure modes and 
effect analysis concept helps to optimize sensor combination based on 
sensor network ability in detecting most probable failures with high 
severity. 

A case study of an air handling unit is presented to illustrate the 
procedure of optimization. Results present optimal sensor combinations 
for different available budgets. Additionally, the application of the 
proposed fault diagnosis method in detecting root causes of abnormal 
energy consumption is briefly discussed. 

According to the results, this study will be helpful for developing 
sensor configuration optimization of any other energy systems. In 
addition, this method could be beneficial in building energy 
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Fig. 15. Scaled information value of failure modes for six optimal sensor combinations. 
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Fig. 16. Actual energy consumption deviation from calculated standard value. 
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Fig. 17. Analyzing root cause of abnormal energy consumption of a signal. 


consumptions by predicting faults in the system. 
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